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WIDEBAND COMMUNICATION 
INTERCEPT AND DIRECTION 
FINDING DEVICE USING HYPERCHANNELIZATION 



The invention disclosed and claimed in this application was developed in 
connection with a contract with the U.S. Government, who as a consequence may have 
rights in the invention. 

Background of the Invention 
I. Field of the Invention 

The present invention relates in general, to communication, and in particular, to 
a communication intercept device using digitally programmable devices to channelize, 
detect, and recognize arbitrary signals. The present invention involves the use of a new 
and improved variable delay First-In-Tap-Out (FITO) spectral memory storage, a two- 
stage synthesis filter for improved dynamic range and programmable software digital 
receivers that use spectral synthesis rather than numerically controlled oscillators for 
signal demodulation. 

In addition, the present invention relates to detection, demodulation, and 
determination of direction of signal activity, and in particular, to a signal intercept and 
direction finding device using general purpose commercial-off-the-shelf (COTS) digitally 
programmable devices to hyperchannelize radro or acoustical spectrum into spectral cells 
narrower than the signal structure of interest, detect, synthesize, and recognize arbitrary 
signal activity, and determine signal direction by phase determination of angle of arrival 
or time difference of arrival. \ 

Overall, the present invention relates to a new and improved multiple and 
cascadable module device using general purpose digital processing, organized in a 
fashion that each module of respective digital signal processing means comprises data 
buffering, digital signal processing and high-speed data routing devices. 
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II. Discussion of the Related Art 

A communication intercept device typically employs a wideband receiver that 
amplifies and down converts a large portion of the radio frequency spectrum. For the 
detection of narrowband signals, the wideband receiver output is channelized by some 
means such as acousto-optic Bragg cells and photo detectors, surface acoustic wave 
(SAW) devices and programmable digital Fast Fourier Transformation (FFT) devices. 

The digital FFT based channelizer is more flexible than the acousto-optic or SAW 
device, and is capable of providing complex (real and imaginary) spectral segments 
sampled above the Nyquist rate. In many applications the channelization segments 
match the bandwidth of the intended receive signal, and each channelization segment is 
therefore connected to a digital signal processing (DSP) device. For example, a series 
of U.S. PATENTS, Nos. 5,535,240, 5,537,435 and 5,848,097 commonly entitled 
"TRANSCEIVER APPARATUS EMPLOYING WIDEBAND FFT CHANNELIZER 
AND INVERSE FFT COMBINER FOR A MULTICHANNEL COMMUNICATION 
NETWORK" disclose a transceiver apparatus that employs this fundamental architecture. 

To reduce the number of DSP devices (e.g. digital drop receivers) required in an 
intercept system, strategies are used to find the spectral channels of activity. For 
example, U.S. PATENT No. 5,428,667, entitled "MULTI-CHANNEL CELLULAR 
COMMUNICATIONS INTERCEPT SYSTEM" discloses a device wherein one or more 
of the plurality of channelized outputs are selected for continuous demodulation. The 
selected channels are the cellular forward control channels that contain formatted 
instructions to mobile radios, establishing the immediate use of narrowband duplex 
channels for cellular network communication. The intent of this invention is to determine 
the cell channel assignments from the control channels and assign only a small plurality 
of drop receivers to the cell channels in immediate use. Estimates of the number of drop 
receivers required to satisfy the demand to a level of confidence can be determined by 
Poisson statistics based on channel activity and cellular call duration as those skilled in 
the art are well aware. 



Further enhancements to wideband intercept systems have been made. U.S. 
PATENT No. 5,812,609, entitled "COMMUNICATION INTERCEPT DEVICE USING 
DIGITAL DROP RECEIVERS IN MULTIPLE TIERS" discloses a device that uses a 
plurality of digital drop receivers in a configuration using a plurality of tiers. This device 
uses fewer digital drop receivers than one per voice grade channel but allows the 
monitoring of all narrowband signals of interest. As exemplified in the device 
description, the first analog to digital converter (ADC) is connected to a first memory that 
stores the digitized wideband signal. A plurality of first digital drop receivers is 
connected to the first memory for selecting various wideband signals in a multiplexed 
fashion under control of a controller. The output of a first digital drop receiver is 
connected to a second memory for storing the medium band signals of interest. The 
second memory has an input connected to the output of each of said first plurality of 
digital drop receivers, storing the signals selected. A second plurality of digital drop 
receivers in a second tier is connected to the second memory for further selection of 
various narrowband signals of interest in a multiplex fashion under control of the 
controller. 

Additional enhancements have been made to wideband communication devices, 
including dynamic signal routing and wideband gain control. U.S. PATENT No. 
5,590,156 entitled "MULTICHANNEL WIDEBAND DIGITAL RECEIVER MAKING 
USE OF MULTIPLE WIDEBAND TUNERS HAVING INDIVIDUALLY 
SELECTABLE GAINS TO EXTEND OVERALL SYSTEM DYNAMIC RANGE" 
discloses a wideband digital receiving device that processes and converts from analog to 
digital representation the bandwidth in a plurality of sub-bands. The output of each sub- 
band is channelized using a complex FFT engine and convolution filtering. According 
to the device description, the convolution filter may be a so-called overlap and add digital 
filter, or may be a polyphase digital filter. Each digital channel signal corresponding to 
a digitized version of the communication signal transmitted by one of the subscriber 
units. Before passing the individual digital channel signals to the DSP processors, each 
digital channel signal output may also be processed by a sample rate adjuster that 



provides samples of its respective digital channel signal taken at or near a position of 
peak symbol amplitude. Other aspects of the invention describe amplifier attenuation 
factors on the wideband receivers based on received signal strength indicator (RSSI) 
made at the output of the channelizer by the DSP processors. As part of the 
communication network, forward channel control reassignment is made to subscriber 
frequencies, grouping like RSSI amplitude signals in frequency channels of a common 
wideband receiver, thereby maximizing channel path ADC performance. 

Multichannel wideband digital receivers are extended to perform direction finding 
by angle of arrival. U.S. PATENT No. 5,815,117 entitled "DIGITAL DIRECTION 
FINDING RECEIVER" discloses a multichannel device capable of receiving radar pulse 
signals simultaneously from a plurality of antennas. Each respective antenna signal is 
received and analog to digital converted. A first antenna digital channelizer divides the 
received antenna signal into N digital frequency channels covering the received spectrum. 
A threshold detector means determines pulse signal presence in a respective channel and 
tunes a digital local oscillator (LO) to the channel frequency as reference. A plurality of 
digital down converters receive respective plurality of delayed antenna digital signals 
from a First in First Out (FIFO) buffer and respective tuning digital LO reference. The 
plurality of digital downconverters output tuned and filtered narrowband signals that are 
phase coherent. The narrowband antenna signal outputs are received by a phase detector 
that compares phase of the first antenna narrowband antenna signal output to all 
subsequent antenna narrowband antenna signal outputs. The plurality of signal phase 
angles and respective signal amplitudes is used in a conventional manner by an angle of 
arrival encoder to estimate the direction of a signal in the selected channel. 

Although wideband receivers and channelizers with drop receivers and 
multichannel processors have been used in the past, it will be readily appreciated that 
improvements are needed to: (1) provide improved signal detection and direction finding 
of signals whose bandwidth do not match a fixed channel bandwidth, (2) provide 
functional isolation and graceful system degradation by re-routing data from one module 
processing device to another, (3) provide flexible ordering of the module processing 



functions, including detection, demodulation and angle of arrival processes, and (4) 
provide efficient digital signal processing using general purpose digital processors 
devices for lower cost and long term supportability. 



Summary of the Invention 

It is an object of the present invention to intercept a wideband RF signal 
environment, sub-band and channelize said signal into a plurality of complex spectral 
segments and select therefrom a plurality of narrowband signals in a flexible manner by 
combining spectral segments in time and frequency for improvement in signal detection 
and demodulation of arbitrary narrowband modulation types and bandwidths. 

It is another object of the present invention to channelize a wideband RF signal 
environment and select therefrom a plurality of narrowband signals using a spectral 
detector employing short-term, medium-term, and delayed long-term spectral averages 
to determine the time of onset of each new narrowband signal. 

It is another object of the present invention to initiate a demodulation process on 
complex spectral data for each new narrowband signal detected using variable delayed 
data from memory storage in a manner that the demodulation process starts at the 
beginning presence of each detected narrowband spectral energy. 

More particularly, in accordance with the present invention, the time-overlapped 
FFT-based channelizers employed in their respective receiver sub-band paths do not 
make use of time-window filters or convolution filters of any kind, thus allowing a choice 
of filtering at later stages. The detection process employs convolution filtering (e.g. 
Nuttal) of the complex spectral data. The demodulation process requires higher dynamic 
range and employs a two-stage synthesis filter (e.g. overlap and add) before complex time 
domain data (i.e. filtered and transformed complex spectral data) is output to a plurality 
of software digital receivers for demodulation. 



It is another object of the present invention to use a plurality of software digital 
receivers for demodulation of certain types of narrowband signals (e.g. two software 
digital receivers may be employed respectively to demodulate the mark and the space of 
a BFSK signal). This demodulation process allows significant elimination of interfering 
signal energy by precise tuning of narrowband filters. 

Brief Description of the Drawings 

Figure la diagrammatically illustrates a single sensor embodiment of the 
invention, channelizing wideband sensor input into a plurality of sub-bands and channels 
for signal processing with sensor antenna 5 connected to a wideband receiver module 
10 comprising wideband receiver 110 connected to a sub-band module 20 comprising 
analog to digital converter 210, a plurality of digital down converters 220, collective 
DMA buffer 230 and output data router port 270 connected to channelizer module 30 
comprising a plurality of input data router ports 307 and respective DMA buffers 308 that 
output to a plurality of FFT channelizers devices 310 and respective output DMA buffers 
330 and output data router ports 370 connected to signal processing module 40 
comprising a plurality of input data router ports 407, a plurality of channel processors 
410 connected to collective data output bus 83 connected to digital recording means 90 
under the control of a system controller 80 controlling apparatus by means of control line 
81 and system bus 82 and data router ports 270, 307, 370, 407 interconnected by means 
of data router device 70 (not shown); 

Figure lb diagrammatically illustrates an N-sensor embodiment of the invention, 
channelizing N wideband sensor inputs into a respective plurality of sub-bands and 
channels for signal processing. Such a variation of the invention is useful for parallel 
comparison of time of arrival of signals for Time Difference of Arrival geolocation 
means; 

Figure 2a diagrammatically illustrates output of time series data from data router 
port 307 and respective DMA buffer 308 to a respective overlap FFT channelizer device 
310 comprising three FFT engines producing a complex spectra output connected to a 



respective output DMA buffer 330 and output high-speed data router port 370 according 
to the present invention; 

Figure 2b schematically illustrates respective three FFT engines FFT-1, FFT-2 
and FFT- 3 of FFT channelizer 310 producing output of time overlapped samples of 
complex spectra; 

Figure 2c schematically illustrates the band filtered RF spectral data, the 
segmented RF complex spectral data from first respective FFT engine, and segmented 
RF complex spectral data from collective FFT-1 to FFT-3 output after filtering; 

Figure 3 diagrammatically illustrates output of complex spectral data segment 
upper half of spectrum band from high-speed data router port 407 and output of complex 
spectral data segment lower half of spectrum band from high-speed data router port 407' 
connected to respective channel processor 410 and 410'. Respective channel processors 
410 and 410' comprise respective FITO delay memory 412, 412', respective detection 
processor 414, 414', respective recognition processor 415, 415' connected to respective 
DMA buffer 430 and 430' for output to digital recorder device 90 by means of data bus 
83 with respective real-time computer controller 413,413' connected to system controller 
80 by means of system bus 82; 

Figure 4 diagrammatically illustrates a respective detection processor 414 and 
respective one demodulation/recognition processor 415 comprising two-stage synthesis 
filter 416 and plurality of synthesis filter receivers 417 under control of synthesis control 
and demodulation processor 435 according to the present invention; 

Figure 5 diagrammatically illustrates a variation of the arrangement of synthesis 
filters in cascading stages according to the present invention wherein first tier of filtering 
is implemented by two-stage filter 416 connected to cascading series of single synthesis 
filters 417; 



Figure6 diagrammatically illustrates N-sensor embodiment of invention to 
perform signal demodulation, recognition, and direction finding. First channel sensor 5 
comprises first channel receiver module 10 and processing chain of respective first 
channel sub-band decimation module 20, channelizer module 30, and signal processor 
module 40. Remaining plurality of sensors comprise a similar plurality of respective 
receiver modules 10'... 10", sub-band decimation modules 20' ...20", and channelizer 
modules 30' ... 30". Collective first channel and remaining sensor channels of 
channelized data is collectively received by DF processing module 60 comprising DF 
processors 610 for respective sub-band decimated channels of each collective N- sensors 
and output data parameters to data bus 83 to recording means 90; 

Figure 7 diagrammatically illustrates alternative embodiment of present invention 
to receive, demodulate, recognize, and DF signals from 2-Channel arrangement. Three 
sensors 5, 5\ 5" representing reference, spatial sine-phase, and spatial cosine phase 
antennas or acoustic sensors are divided into to processing chains. First channel sensor 
5 comprises first channel receiver of dual channel receiver 111 of module 10 and 
processing chain of respective first channel sub-band decimation module 20, channelizer 
module 30, and signal processor module 40. Sensor output 5, 5', 5" are cyclically 
sampled by commutation switch and output to second receiver 1 1 1 of receiver module 
10 connected sub-band decimation module 20', and channelizer module 30'. Collective 
first channel and commutated second channel channelized data is collectively received 
by DF processing module 60 comprising DF processors 610 for respective sub-band 
decimated channels and output data parameters to data bus 83 to recording means 90; 

Figure 8a illustrates a variation of the DF processor of present invention for 
receiving N-sensor hyperchannel spectral input. DF processor comprises multiplexer 
611, 61 1 5 for receiving upper and lower spectral half data and connected to respective 
First In Tap Out (FITO) delay memory buffers for collective output to respective DF 
synthesis filters 614, 614', that filter and transform hyperchannel spectral data for 
collective output to respective DF calculators 616, 616' and angle of arrival output to 
system data bus 83; 
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Figure 8b illustrates a variation of the DF processor of present invention for 
receiving 3-Channel reference, spatial sine-phase, and spatial cosine-phase sensor data. 
DF processor comprises multiplexer 61 1, 611 ' for receiving upper and lower spectral half 
data and connected to respective First In Tap Out (FITO) delay memory buffers for 
collective output to respective DF synthesis filters 614, 614', filter and transform 
hyperchannel spectral data for collective output to respective DF calculators 616, 616' 
and angle of arrival output to system Data Bus S3. 3 -Channel signal demodulation, 
recognition, and DF is an embodiment of the N-sensor embodiment of present invention 
of particular usefulness; 

Figure 8c illustrates a variation of the DF processor of present invention for 
receiving 2-Channel hyperchannel spectral input. DF processor comprises multiplex 611, 
611' for receiving upper and lower spectral half data, performing decommutation and 
calibration of second channel commutated data input, forming 3 -Channel plieso- 
synchronous reference, spatial sine-phase, and spatial cosine-phase hyperchannels for 
output to respective three upper and three lower spectral half FITO delay memory buffers 
612, 612' for collective output to respective upper and lower half DF synthesis filters 
614, 614% that filter and transform three-channel hyperchannel spectral data to time 
domain data for collective output to respective DF calculators 616, 616' and angle of 
arrival output to system data bus 83; 

Figure 8d illustrates a variation of the present invention for beam forming INT- 
sensor hyperchannel spectral input. Beamforming processor 690 comprises multiplexer 
61 1, 61 1' for receiving upper and lower spectral half data connected to respective First 
In Tap Out (FITO) delay memory buffers for collective output to respective N-sensor 
beamformers 617, 617' that filter apply beamforming weights to each N-sensor 
hyperchannel input, providing one or more plurality weighted beamformed outputs to 
respective high speed data port 670, 670'; 

Figure 9a illustrates an N-sensor two-stage DF synthesis filter 414 comprising 
first stage multipliers 621a to 62 In receiving hyperchannel data from respective FITO 



buffers 612a to 61 2n and multiplying input by common filter weight 620 with common 
shifted indices set by step tune selector 622 output to respective first stage inverse FFT 
devices 623a to 623n. Inverse FFT devices are 623a to 623n are connected to respective 
contra-phase I&Q rotation devices 624a to 624n with output to respective first stage 
FIFO buffers 625a to 625n. Second stage respective FFT devices 626a to 626n receive 
output from respective first stage FIFO buffers 625a to 625n and in manner similar to 
first stage are filtered by respective multipliers 627a to 627n with common filter weights 
620\ Respective hyperchannel products are inverse transformed by respective FFT 
devices 623'a to 623 'n connected to respective FIFO buffers 625'a to 625'n for output; 

Figure 9b illustrates an N-sensor beamformer 617 comprising first stage 
multipliers 621a to 62 In receiving hyperchannel data from respective FITO buffers 612a 
to 612n and multiplying input by common filter weight 620 with common shifted indices 
set by step tune selector 622 output to respective first stage inverse FFT devices 623a to 
623n. Inverse FFT devices are 623a to 623n are connected to respective contra-phase 
I&Q rotation devices 624a to 624n with output to respective first stage FIFO buffers 625a 
to 625n. Tuned and filtered hyperchannel time-domain output of FIFO buffers 625a to 
625n is received by respective multiplier devices 631a to 63 In that multiply respective 
N-sensor time domain data by N-sensor complex weights 630 connected to summer 
device 635 that collectively adds weighted N-sensor data into a single beamformed time 
domain output to FIFO buffer 636. Time domain beamformed output 636 is transformed 
to frequency domain by FFT 637 for output to high speed data router 70. It is to be 
understood that time domain data may be selected tor output in additional variations of 
N-sensor beamformer 617; 

Figure 1 0a illustrates data flow through N-sensor multiplexer 61 1 of DF processor 

610 and Beamform processor 690. Respective hyperchannel data from each N-sensor 
FFT hyperchannelizer (upper or lower spectral half) is received by respective multiplexer 

6 1 1 by means of high speed data router and ports (not shown). Multiplexer orders 
respective hyperchannel data from respective N-sensors to respective FITO delay 
memory buffer; 
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Figure 10b illustrates data flow through 3-channel multiplexer 611 of DF 
processor 610. Respective hyperchannel data from reference, sine-phase, and cosine 
phase hyperchannel FFT hyperchannelizer (upper or lower spectral half) is received by 
respective multiplexer 61 1 by means of high speed data router and ports (not shown). 
Multiplexer orders respective hyperchannel data from respective reference, sine-phase, 
and cosine phase to respective FITO delay memory buffer; 

Figure 10c illustrates data flow through 2-channel multiplexer 611 of DF 
processor 610. Respective hyperchannel data from first channel reference hyperchannel 
data and second channel commutated reference, sine-phase, and cosine-phase 
hyperchannel data from respective FFT hyperchannelizer (upper or lower spectral half) 
is received by respective multiplexer 61 1 by means of high speed data router and ports 
(not shown). Multiplexer orders respective hyperchannel data of first channel reference 
data and second channel commutated data. Multiplexer compares time related reference 
data from first and second channel for calibration. Difference between respective 
reference value between first and second channel is added to respective second channel 
reference, sine-phase and cosine-phase data of second channel. Parallel output of second 
channel reference, sine-phase and cosine-phase data is at one third the input data rate; 

Figure 11a illustrates one variation of the ordering of processing modules of the 
present invention for a single sensor system to perform signal activity detection and 
recognition, Sensor 5 connects to receiver module 10 providing analog wideband IF 
output to sub-band module 20 that performs analog to digital conversion and sub-band 
decimation connected to channelizer module 30 by high speed router 70. Channelizer 
module 30 performs respective sub-band FFT and outputs respective hyperchannelized 
spectral data to signal processor module 40 connected by high speed router 70. Signal 
processor module 40 performs signal activity detection and spawns signal demodulation 
and recognition processors using synthesis filtering and parallel demodulators to 
determine signal of interest. Signal of interest parameters are output by means of data 
bus 83 to digital recorder means 90. System controller 80 controls receivers by means 
of control line 8 1 and digital modules by means of system bus 82; 
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Figure 1 lb illustrates another variation of the ordering of processing modules of 
the present invention for an N-sensor system to perform signal activity detection, 
recognition, and direction finding. Sensors 5 to 5" connect to respective receiver 
modules 10 to 10" providing analog wideband IF output to respective sub-band modules 
20 to 20" that perform analog to digital conversion and sub-band decimation connected 
to channelizer modules 30 to 30" by high speed router 70. Channelizer modules 30 to 
30" perform respective sub-band FFT and outputs respective hyperchannelized spectral 
data to respective signal processor modules 40 to 40" connected by high speed router 70. 
First channel signal processor modules 40 performs signal activity detection and spawns 
signal demodulation and recognition processors using synthesis filtering and parallel 
demodulators to determine signal of interest. Signal of interest parameters are output to 
DF processor 60 that receives collective N-sensor hyperchannel data from FFT 
channelizer modules 30 to 30" and outputs angle of arrival parameters for respective 
signal of interest to data bus 83 connected to digital recording means 90. System 
controller 80 controls receivers by means of control line 81 and digital modules by means 
of system bus 82; 

Figure 1 lc illustrates another variation of the ordering of processing modules of 
the present invention for a 3-sensor system commutated system to perform signal activity 
detection, recognition, and commutated direction finding. Reference sensor 5 connects 
to respective first receiver of dual receiver 1 1 1 providing analog wideband IF output to 
respective first sub-band module 20. Reference sensor 5, sine-phase sensor 5 5 and 
cosine-phase sensor 5" connect to commutative switch 112 that cyclically outputs 
reference, sine-phase, and cosine-phase RF or acoustic signal to second receiver of dual 
receiver 1 1 1 providing analog wideband IF output to respective second sub-band module 
20" that perform analog to digital conversion and sub-band decimation. Respective first 
and second sub-band modules 20 and 20" are connected to respective channelizer 
modules 30 and 30" by high speed router 70. Respective channelizer module 30 and 30" 
perform respective sub-band FFT and outputs respective hyperchannelized spectral data 
to respective signal processor modules 40 and 40" connected by high speed router 70. 
First channel signal processor modules 40 performs signal activity detection and spawns 
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signal demodulation and recognition processors using synthesis filtering and parallel 
demodulators to determine signal of interest. Signal of interest parameters are output to 
DF processor 60 that receives first channel reference and second channel cyclic reference, 
sine-phase, and cosine-phase hyperchannel data from respective FFT channeiizer 
modules 30 and 30". DF processor 60 comprising multiplexer 611 (not shown) performs 
calibration of second channel reference, sine-phase, and cosine phase data and outputs 
hyperchannel data in parallel through synthesis filters (not shown) and collective DF 
calculator (not shown) and outputs angle of arrival parameters for respective signal of 
interest to data bus 83 connected to digital recording means 90. System controller 80 
controls receivers by means of control line 81 and digital modules by means of system 
bus 82; 

Figure 1 Id illustrates another variation of the ordering of processing modules of 
the present invention for an N-sensor system to perform direction finding, constant angle 
of arrival activity, and signal of interest recognition. Sensors 5 to 5" connect to 
respective receiver modules 10 to 10" providing analog wideband IF output to respective 
sub-band modules 20 to 20" that perform analog to digital conversion and sub-band 
decimation connected to channeiizer modules 30 to 30" by high speed router 70. 
Channeiizer modules 30 to 30" perform respective sub-band FFT and outputs respective 
hyperchannelized spectral data to respective collective DF processor module 60 and to 
first channel signal processor module 40 connected by high speed router 70. DF 
processor module performs N-sensor direction of arrival calculation on each 
hyperchannel with angle of arrival parameters output to first channel signal processor 
module 40 by means of high speed data router 70. Detection processor 414 (not shown) 
of first channel signal processor module performs direction of arrival sorting, and spawns 
signal demodulation and recognition processors 415 (not shown) using synthesis filtering 
to tune to each constant angle hyperchannel and using parallel demodulators to determine 
signal of interest. Signal of interest parameters are output to data bus 83 connected to 
digital recording means 90. System controller 80 controls receivers by means of control 
line 81 and digital modules by means of system bus 82; and 
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Figure lie illustrates another variation of the present invention for an N-sensor 
system to perform beamforming signal to noise enhancement, signal activity detection, 
recognition, and direction finding. Sensors 5 to 5" connect to respective receiver 
modules 10 to 10" providing analog wideband IF output to respective sub-band modules 
20 to 20" that perform analog to digital conversion and sub-band decimation connected 
to channelizer modules 30 to 30" by high speed router 70. Channelizer modules 30 to 
30" perform respective sub-band FFT and outputs respective hyperchannelized spectral 
data to collective beamformer module 690 and to collective DF processor module 60 by 
means of high speed data router 70. Beamformer 690 performs one or more plurality N- 
sensor weighting and N-sensor hyperchannel combining, outputting one or more plurality 
of spatially directive beams of hyperchannel data to respective plurality of signal 
processor modules 40 by high speed router 70. Signal of interest parameters are output 
to respective beam DF processor 60 that receives collective N-sensor beam hyperchannel 
data from channelizer module 30 and output angle of arrival parameters for respective 
signal of interest to data bus 83 connected to digital recording means 90. System 
controller 80 controls receivers by means of control line 81 and digital modules by means 
of system bus 82. 

Detailed Description of the Preferred Embodiments 

With reference to the figures, like reference characters will be used to indicate like 
elements throughout the several embodiments and views thereof. Before describing in 
detail the particular improved wideband multichannel receiving device in accordance 
with the present invention, it should be observed that the present invention resides 
primarily in a novel combination of commercially available communication and signal 
processing circuits and components and not in the particular detailed configurations 
thereof. Accordingly, the structure, control, and arrangement of these conventional 
circuits and components have been illustrated in the drawings by readily understandable 
block diagrams that show only those specific details that are pertinent to the present 
invention, so as not to obscure the disclosure with structural details which will be readily 
apparent to those skilled in the art having the benefit of the description herein. Thus, the 
block diagram illustrations of the Figures do not necessarily represent the mechanical 
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structural arrangement of the exemplary system, but are primarily intended to illustrate 
the major components of the system in a convenient functional grouping, whereby the 
present invention may be more readily understood. 

In addition, as one of skill in the art will understand from this disclosure of the 
5 invention, some components of the system may be implemented by hardware, other 

components may be implemented by software, while still other components may be 
implemented by either software or hardware, or even a combination of both. 
Specifically, as will be explained in more detail hereinafter, the invention is composed 
of elements that will occur or exist in pluralities (i.e., digital down converters, FFT 

10 channelizers, channel processors, FFT engines, recognition processors). One of skill in 

the art would understand that, if implemented as hardware, the embodiment of such plural 
components would involve parallel sets of such components that would be selectively 
activated depending upon the specific application of the system. If implemented as 
software, the embodiment of those components would be the activation of a subroutine, 

15 applet, etc. whereby the specific application of the system would determine the number 

of times that the subroutine is repeated or run simultaneously. In any event, such 
variations to the implementations of the invention would all be within the scope of this 
application and the claims therein. 

Wideband Receiver 

20 Referring now to Figure la, the apparatus of the present invention is 

diagrammatically illustrated receiving radio frequency signals by means of an antenna 
5 and wideband receiver 1 10. For convenience, the antenna source, and analog filtering 
circuitry are not shown. The wideband receiver 110 amplifies and down-converts to 
baseband a time domain analog signal containing the contents of all radio frequency 

25 signals received in a selected wideband band-pass. As a non-limiting example, the 

wideband receiver 1 10 may comprise a a high dynamic range receiver model 9 104 A, 
manufactured by Watkins Johnson with analog output of 1 0 MHZ in the High Frequency 
(HF) Spectrum (1-30 MHZ). Another example of a wideband receiver is the VXI-3570 
manufactured by Cubic Cubic Corporation with tunable analog output of 4 MHZ 
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bandwidth in the Very High Frequency (VHF) Spectrum (30-300 MHZ) and Ultra High 
Frequency (UHF) Spectrum (300-3000 MHZ). The wideband receiver 1 10 is tuned from 
a local control line 81 connected to the system controller 80. 

The wideband radio frequency baseband output from the wideband receiver 1 1 0 
is coupled to the high-speed analog to digital converter (ADC) 210 in the sub-band 
module 20. Advantageously, the dynamic range and sampling rate capability of currently 
commercially available ADCs, such as that referenced above, are sufficiently high that 
the signal environment is linearly represented over 80dB dynamic range. 

It should be observed that the present invention is not limited to 10 MHZ 
bandwidth or any other set of system parameters. The values given here are merely for 
the purposes of providing an illustrative example. Also, while the term 'wideband 5 is not 
limited to any particular spectral range, it is to be understood to imply a spectral coverage 
of at least an order of magnitude larger that the bandwidth of the majority of signals in 
that portion of the spectrum. Narrowband, on the other hand, implies only a portion of 
the spectrum, for example, the width of an individual channel (e.g. 3 kHz in HF, 30 kHz 
in VHF or UHF). 

Sub-Band Module 

The output of wideband receiver 1 10 of the receiver module 10 is received by the 
Sub-Band Module 20, which incorporates the analog-to-digital converter (ADC) 210, a 
plurality of Digital Down Converters (DDCs) , a collective demand memory access 
(DMA) digital data buffer 230 and respective high-speed data router port 270. 

Typical operation of an ADC would be slightly above the Nyquist sampling rate 
to ensure that the signal data is not aliased. For the present invention, the ADC 210 is 
operated significantly higher than the Nyquist sampling rate (e.g. for a wideband receiver 
bandwidth of 10 MHZ, the ADC sampling rate is 32.4 MHZ), ensuring that the complex 
signal is oversampled in time. This is an important feature of the present invention, 
allowing subsequent channelization, and channel recombination in time and frequency 
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to occur without distortion. Another important feature of the present invention, allowing 
subsequent time of arrival measurements, is time synchronization of ADC operation, 
allowing data to be '"time-tagged". The digital complex in-phase and quadrature (I&Q) 
oversampled signal data output is distributed to the plurality of Digital Down Converters 
5 (DDCs) 220 for sub-band tuning. In one embodiment, the DDCs 220 may be directly 

connected to the ADC 210 in a mezzanine assembly of a card chassis assembly (CCA). 
Alternatively, the time domain digital data from the ADC 210 may be received by a data 
buffer (not shown) and then outputted asynchronously to the DDCs 220. 

In general, the DDCs 220 receive the analog-to-digital wideband data from the 
10 ADC 210 and output at or near baseband a respective small plurality of sub-band time 

domain complex digital data output for distribution by the high-speed data router port 

1*1 

*0 270. In the preferred embodiment of the present invention, for a wideband receiver 110 

ill with 1 0 MHZ bandwidth, the digital radio frequency signal is distributed to six DDCs 

y t 220, each respective ones DDC providing a 1.8 MHZ sub-band complex time-domain 

W is signal. The sub-band tuning of the DDCs 220 is determined by the system controller 80 

■P 

connected to the respective ones DDCs by a system control bus 82 (e.g. VMEbus™). 

j"* The digitized data complex time domain signals output by the plurality of DDCs 

if! 220 are collectively buffered in a demand access memory (DMA) buffer 230 for access 

by a high-speed data router port 270 using the principles of Time Division Multiplex 
20 (TDM) access to the memory of the collective buffer 163. The collective buffer 163 may 

be implemented, for example by Application Specific Integrated Circuit (ASIC) 
components supplied on the EchoTek digital down converter card, Model ECDR-412. 

High-speed TDM access of the data router port 270 to the memory data of the 
collective buffer 230 with data routing at clock rates substantially higher than the 
25 collective DDC data rates are output to one or more of a plurality of destination high 

speed data router ports. The TDM access and the data routing are configured by the 
system controller 80 that is connected to the high-speed data router for control of the 
respective data router ports by the system control bus 82 (e.g. VMEbus™). One-to-one 
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or one- to-many high-speed data links may be implemented using, for example, Race Way 
Interlink model ILK6 manufactured by Mercury Computer Systems, Inc. 

In the preferred embodiment of the invention, high-speed data routing links are 
configured in a one-to-one arrangement from the plurality of DDCs 220 and respective 
5 collective DMA buffer 230 to a respective plurality of high-speed data router ports 307 

for output to a plurality of FFT channelizers 310 in the Channelizer Module 30. 

Again, it is to be understood that the term 'sub-band' is not limited to any 
particular spectral range, but imply a signal bandwidth larger that the bandwidth of the 
majority of signals in that portion of the RF or acoustic spectrum. In one embodiment 
10 of the invention, for HF radio data as an example, the digital down converters 220 would 

be configured to output six decimation sub-bands. Each respective sub-band would 
comprise near baseband complex time-series data at 2.7 mega-samples per second. 

High-Speed Data Router 

Time Division Multiplex (TDM) memory access is performed by high-speed data 

15 router ports to data at the collective buffer 230 and data routing are at clock rates 

substantially higher than the collective DDC 220 data rates. Data routing is configured 
by the system controller 80 that is connected to the high-speed data router 70 (not shown) 
and respective port 270 by the system control bus 82 (e.g. VMEbus™). It is to be 
understood that the high-speed data linkage is implemented by means of non-blocking 

20 data router ports that digitally interconnect modules according to various topologies. 

In the preferred embodiment of the invention, high-speed data routing of one-to- 
one and one-to-many may be configured. Referring to Figure la, the digitized data 
complex time domain signals of first DDC sub-band tuner 220 of sub-band module 20 
is buffered by collective buffer 230 and routed from high-speed data port 270 by Time 
25 Division Multiplexing (TDM) to first high-speed data port 307 for asynchronous data 

buffering by first demand memory access (DMA) buffer 308 and output to first FFT 
channelizer 310 of channelizer module 30. 
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Further, non-blocking routing of one-to-many is illustrated. Referring to Figure 
la, digital data complex time domain signals from respective plurality of DDC sub-band 
tuner components 220 are buffered by collective buffer 230 and routed from high-speed 
data port 270 by Time Division Multiplexing (TDM) means to respective plurality high- 
5 speed data ports 307 for asynchronous data buffering by respective plurality of DMA 

buffers 308 and output to respective plurality of FFT channelizers 310. 

In the embodiment of the preferred invention, high-speed data routing links may 
be configured in a one-to-one or one-to-many arrangement from a P plurality of FFT 
channelizers to a Q plurality of channel processors. When configured for one-to-one 
10 links between respective FFT channelizer 310 and output to respective DMA buffer 330 

for transfer by respective high-speed data router port 370, respective TDM data transfer 
is made to respective high-speed data router port 407 and output to respective channel 
processor 410. 

When configured, for example, for one-to-a plurality of links between one FFT 
1 5 channelizer and a plurality of channel processors, in the embodiment of the preferred 

invention, two channel processors are employed. Referring to Figure 3, in such an 
arrangement, upper and lower halves of complex spectral data segments are output by an 
exemplified FFT channelizer 310' and DMA buffered 330' for asynchronous TDM 
transmission by respective high-speed data port 370'. High-speed data router port 370' 
20 separates complex spectral data segments into upper and lower complex spectral data 

segments, and transmits respective upper half and lower half spectral segments on links 
TDM-1 and TDM-2. Link TDM-1 of high-speed data router port 370' outputs upper half 
spectral data segments to respective high-speed data router port 407' providing 
asynchronous upper half spectral data segments to channel processor 410'. Link TDM-2 
25 of high-speed data router port 370 outputs lower half spectral data segments to respective 

high-speed data router port 407" providing asynchronous lower half spectral data 
segments to channel processor 410". 
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Non-blocking routing of TDM high-speed data in one-to-one or one-to-many 
transfer links may be implemented using, for example, Race Way Interlink data router 
Model ILK6 manufactured by Mercury Computer Systems, Inc. 

FFT Channelizer 

5 The channelizer module 20 is designed to receive time domain complex digital 

data from the high-speed data router port 270 of the sub-band module 20 to input high- 
speed data router port 307 connected to a DMA data buffer 308 of a respective FFT 
hyperchannelizer 310. The FIFO data buffer 308 performs buffering for asynchronous 
data arrival and outputs complex time domain data to its corresponding respective Fast 
1 0 Fourier Transform (FFT) device 310. 

*0 The FFT device 310 comprises one or more plurality of FFT engines (FFT1 - 

ill FFT3) producing complex spectra that are overlapped in time. In the present 

embodiment, the FFT device 310 may select complex transforms 2 A N in size. Unlike 
conventional FFT processes, the FFT device 310 is configured to perform no spectral 
»! 15 windowing. For illustration, an 8k FFT results in a spectral data bandwidth of 

\\l approximately 330 Hz that is 10% the bandwidth of expected signals. The object of the 

!** FFT device 3 1 0 is to channelize time domain data into complex spectral components that 

*£! are much narrower than any signal of interest bandwidth. This is referred to as 

hyperchannelization. The output of the FFT device 310 is connected to a respective 
20 DMA buffer 330 for asynchronous output to high-speed data router port 370 for 

distribution to a respective spectral detection device or other signal processing module. 
In another embodiment of the invention, if a multichannel antenna or multichannel source 
data is present, the output of the FFT device 310 buffered by DMA buffer 330 may also 
be routed via the high-speed data router port 370 to a direction finding module. In still 
25 further embodiments of the invention, the output of the FFT device 310 buffered by 

DMA buffer 330 may be routed via the high-speed data router port 370 to any other 
signal processing module with respective high-speed data router port input means. 
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Referring to Figure 2a and Figure2b each FFT channelizer 310 comprises an input 
high-speed data port 307 and DMA buffer 308 connected to the FFT device 310 
composed of three parallel FFT engines FFT1 - FFT3, with complex spectral output 
collected by DMA buffer 330 and distributed by a high-speed data router port 370. The 
high-speed data router port 370 is connected to a respective high-speed data router port 
407' of first channel processor 410' of signal processing module 40, and if configured 
by the high-speed data router, to a second high-speed data router port channel processor 
407". If the data router 370 of the channelizer module 30 configures both first and 
second channel processors, the first channel processor 410' receives the upper half of 
FFT channelizer complex spectra segments outputted from a corresponding FFT 
channelizer 310, while the second channel processor 410" receives the lower half of the 
FFT channelizer complex spectra segments. 

Referring to Figure 3, representative channel processor 410' incorporates a real- 
time controller 413', a complex spectra variable delay First-In-Tap-Out (FITO) memory 
412, a respective detection processor 414, and a respective plurality of demodulation 
processors 415 connected to respective DMA buffer 430' and high-speed data router port 
470' for complex signal data output. 

Each respective FFT unit 310 performs an N-point complex FFT transformation 
using a block of complex signal data contained in the DMA buffer 308. Referring to 
Figure2b, the operation of the three FFT units comprising FFT unit 3 1 0 is orchestrated 
in time such that there is an overlap in data used by the first and second FFT units, the 
second and third FFT units, and in a cyclic manner, the third and first FFT units. In the 
embodiment of the preferred invention, the overlap is 12.5 percent. However, the device 
is not limited to such specific overlap values, and other overlaps such as 25 percent may 
be used. 

The cyclic output of the three FFT engines FFT-1, FFT-2, FFT-3 that comprise 
FFT unit 3 10 is collectively buffered in a demand access memory (DMA) buffer 330 for 
asynchronous access by a high-speed data router port 370. 
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In the embodiment of the preferred invention, high-speed data routing links may 
be configured in a one-to-one or one-to-many arrangement from an N plurality of FFT 
channelizers to a P plurality of channel processors. When configured for one-to-one links 
between respective FFT channelizers 310 and channel processors 410, all FFT 

5 channelizer complex channel spectral segment outputs are transferred by the high-speed 

data router port 370 to respective channel processors high-speed data router port 407. 
When configured, for example, for one-to-a plurality of links between one FFT 
channelizer and a plurality of channel processors, in the embodiment of the preferred 
invention, two channel processors are employed. In such an arrangement, the respective 

10 first channel processor 410' receives by means of the high-speed data router port 407' the 

complex spectral segments output from the upper half of the respective FFT channelizer 
310 and the respective second channel processor 410" receives by means of the high- 
speed data router port 407" the complex spectral segments output from the lower half of 
the respective FFT channelizer 310. 

15 Channel Processor 

Referring to FIGS. 3 and 4, each respective channel processor 410 incorporates 
a variable First-In-Tap-Out (FITO) delay memory 412 to receive complex spectral data 
from input high-speed data router port 407. FITO delay memory is connected to a 
detection processor 414 and to a plurality of demodulation/recognition processors 415 

20 under control of a real-time controller 413 connected to the system controller 80. The 

detection processor 414 is composed of a shape filter 440, a convolution combiner 441, 
a spectral segment squaring component 442, a spectral parameter encoder 443 a short 
term average buffer AVI , a medium term average buffer AV2, and a delayed long term 
average buffer AV3. 

25 In general, the detection processor 414 conducts spectral convolution of a signal 

inputted from the delay memory 412. The spectral parameter encoder 443 of the 
detection processor, in the presence of new signal energy within the resolution bandwidth 
of the FITO spectral segments, compares the three running averages stored in the buffers 
AVI - AV3, determines the onset of specified signal types as detected by the shape of 
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a plurality of spectral segments and outputs a parameter word containing signal type, 
frequency, time, and bandwidth to the real-time controller 413. 

The real-time controller activates (that is, spawns) the signal 
demodulation/recognition processors 41 5 in response to each new signal parameter word. 
5 Each demodulation/recognition processor 415 is dynamically connected to the FITO 

delay memory 412 to receive the complex spectral data at the start of signal presence. 

The demodulation/recognition processor 415 includes at least a two-stage 
synthesis filter operation employing a first stage synthesis filter comprising band filter 
420 and multiplier 421 circuit with a first offset step tune selector 422, a first inverse FFT 

10 (IFFT) device 423, a first In-Phase and Quadrature (I&Q) signal component rotation 

compensator 424 for tuning, and a first FIFO DMA buffer 425 for asynchronous output 
to second stage synthesis filter. First stage synthesis filter is connected to a second stage 
synthesis filter comprising FFT device 426, a second filter 420', a second multiplier 42 V 
providing complex spectral data output to second IFFT device 423' and a second FIFO 

15 DMA buffer 425' for asynchronous output. 

A second plurality of synthesis filters 416' may be employed for further spectral 
refinement of the signal data. The output of the second FIFO DMA buffer 425' may be 
configured so as to be available for system storage (As illustrated, 425' is configured for 
data output to system storage via high speed data router pout 470). Refined complex time 
20 domain signal data output from the plurality of second FIFO DMA buffers 425' is 

connected to the synthesis receiver controller 434. Synthesis receiver controller 434 
applies demodulated algorithms to collective refined complex time domain signal data 
and output is passed to the real-time controller 413 which passes output to system 
controller 80. 

25 Specifically, the FITO delay memory 412 receives respective FFT channelizer 

complex spectral segment outputs by means of the high-speed data router 407. The FITO 
memory 412 acts in principle as a First-In-First-Out (FIFO) memory storage or delay 
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device, with constant refresh of new data (first in) and constant drop of old data (first out) 
with respect to a determined memory storage allocation size. Unlike typical FIFO 
memory or delay devices, the FITO may be accessed at any relative memory point from 
the first data (that is, tapped out). This provides a variable delay storage access device, 
5 with delayed output ranging from no delay (accessed at first in), to incremental delay 

(accessed at memory points with delay increments equal to the data sampling), to 
maximum delay (accessed at first out, the last point in memory before data is dropped). 
The implementation of the FITO delay memory 412 can be constructed in linear or 
circular data structures as those skilled in the art of data programming would be able to 
10 discern. 

The detection processor 414 accesses respective FFT channelizer complex 
spectral data at the first in (no delay) output point of the FITO memory 412. The 
detection processor 414 comprises the operations of spectral convolution by the shaping 
filter H D (f) 440 and convolution combiner 441, and outputs complex spectra to the 

15 spectral magnitude converter 442 that outputs the spectral magnitude component to the 

parameter encoding 443. Parameter encoder 443 creates, maintains and evaluates short- 
term, medium-term, and delayed long-term spectral averages for each respective real 
spectral segment data. The object design of the shaping filter H 0 (f) 440 is to use a 
convolution filter comprising only real values, thereby reducing the number of multiply 

20 and adds in the convolution process by a factor of two over a truly complex shaping filter. 

More specifically, the convolution processes conducted by the shaping filter 440 
and convolution combiner 441 embodies the process of spectral windowing, wherein 
complex spectral data is convolved with the Discrete Fourier Transform (DFT) of an N- 
term time-domain window function as exemplified, but not limited to, such window 
25 functions described by Gnathal, "Some Windows with Very Good Sidelobe Behavior", 

IEEE Transactions on Acoustics. Speech, and Signal Processing , Vol ASSP-29, No 1, 
Feb 1981. The number of non-zero DFT coefficients is 2N- 1 . In at least one embodiment 
of the present invention, N=4, creating a window function with maximum sidelobe level 
of -93 dB and a hyperchannel bandwidth of approximately twice that of the original 
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respective hyperchannel. The spectral windowing process repeats convolution on 
respective time-decimated sequence of complex spectral data, producing a time-sequence 
of windowed hyperchannels. 



5 The spectral magnitudes generated by the spectral magnitude squaring component 

442 are maintained in running averages of three different lengths, designated the short- 
term, medium-term, and delayed long-term spectral averages for each respective real 
spectral segment in the buffers AVI - AV3. The determination of the lengths of time 
designating short, medium, and long are values established for the particular spectral 

10 environment, with values ranging from, but not limited to, milliseconds to seconds. The 

delayed long-term average is to be understood as an operation that employs a FIFO buffer 
before adding the spectral segment content to the running average. 



The spectral parameter encoder 443 employs the three spectral averages for 
determining the presence of new signal activity and outputs a parameter word including 

15 frequency, time of start, bandwidth, and possible signal modulation type to the real-time 

channel processor controller 413 for activation of one or more demodulation/recognition 
processors 415. It is well known to skilled workers in the art of communication that 
various signal modulations have distinguishing spectral shapes (e.g. distinct peaks with 
spectral separation for Frequency Shift Key (FSK) modulation, and sine(x)/x spectral 

20 humps for Phase Shift Key (PSK) modulation, etc.). 



In essence, in at least one embodiment of the invention, spectral averaging as 
conducted by the combination of the magnitude squaring component 442 and the spectral 
parameter encoder 443 involves receiving time-sequence windowed hyperchannels, 
extracting the magnitude component of each respective complex hyperchannel spectral 
25 value and passes respective magnitude to three simultaneous finite impulse response 

(FIR) lowpass filter and decimation devices. In this case, the FIR devices are embodied 
by the buffers AVI - AV3. As noted above, these FIR devices produce the respective 
short-term, medium-term, and delayed long-term moving averages (STA, MTA, and 
DLTA) of spectra magnitude. In a preferred embodiment, the spectral average data 
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from the buffers AVI - AV3 are then output to high-speed data routing device (not 
shown) to, for example, a visual display device for a real-time spectral display. 

In an alternative embodiment, infinite impulse response (IIR) filters may be used 
to generate the STA, MTA, and DLTA of the buffers AVI - AV3. A typical one-pole 
5 filter computes the Nth output value as a weighted average of the Nth input and the N-l 

average using the recursive algorithm: 

Average (N) = a x Spectra (N) + (1-a) x Average (N-l), where a < 1 

The constant a controls the attack and decay characteristics of the filter. Respective 
alphas are selected for respective short term, medium term, and long term IIR filters. 

1 0 The plurality of demodulation/recognition processors 415 each access a respective 

one of the FFT channelizer complex spectral data at output points in the FITO memory 
412 receiving a delayed data stream in such a manner that the complex spectral data 
received from the FITO memory 412 starts at the beginning presence of a detected signal 
energy. The demodulation/recognition processors 415 each comprise the operations of 

15 synthesis filtering by a two-stage overlap and add filter using only complex multiplies, 

and a plurality of software digital receivers in one or more cascaded stages under control 
of a real time demodulation processor. 

Specifically, the synthesis filtering is accomplished by receiving complex spectral 
segment data from the output of the FITO memory 412 and providing two stages of 

20 filtering using the first spectral band filter Hl(f) 420 employing the first offset multiplier 

42 1 of spectral segment as a means of filter tuning 422 under control of the synthesis 
receiver controller and demodulation processor 435, the first inverse FFT (IFFT) 423 for 
returning the complex spectra to the complex time-domain followed by a first phase 
rotation (I&Q rotate) via the signal component rotation compensator 424 to restore phase 

25 distortion incurred in the first filter and first IFFT transformation due to received FFT 

time-overlapped channelizer data and first filter offset tuning, and the FIFO DMA buffer 
425 for output distribution. 
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The single stage synthesis filter, commonly known to those skilled in the art of 
digital signal processing as an overlap and add filter, is exemplified in Theory and 
Application of Digital Signal Processing . L.R.Rabiner & B.Gold, Prentice-Hall, 1975, 
pp 65-66. The embodiment of the present invention is unique, however, in that a second 
5 stage is subsequently employed for increased dynamic range. The complex time domain 

output of the first synthesis filter stage is connected to a second stage, employing the 
operations the second FFT device 426 (the first being the FFT based channelizer 
producing the time overlapped complex spectra output to the channel processor and 
distributed to the demodulation processor by means of the variable delay FITO memory 

10 412), the second multiplying spectral filter H2(f) 420' by means of second multiply 

device 42 T, Output from second stage multiply device 42 1' is received by second stage 
IFFT device 423' and connected to second stage FIFO DMA buffer 425' for output 
distribution to a plurality of synthesis filters receivers 417. Two stage synthesis filter 416 
and plurality of synthesis filter receivers 417 are under parametric control of the synthesis 

15 receiver control and demodulation processor 435 that receives direction from and passes 

parameters to the real time computer controller 41 3. Additional stages of refinement may 
be employed in the present invention. Referring to FigureS, a preferred embodiment of 
the invention employs additional synthesis filter modules in a cascaded fashion. 



The plurality of software synthesis filters 417 connect to the synthesis receiver 
20 control and demodulation processor 435 to perform signal demodulation. In the 

embodiment of the preferred invention, three or more software digital synthesis filters 
may be simultaneously employed. It is to be noted that for most communication 
purposes, the bandwidth of the receiving digital filter matches the bandwidth of the 
wideband FFT channelizer. Whereas for this present invention, the software digital 
25 filters 417 may act upon a signal that is either wider than or narrower than a single 

respective wideband FFT channelizer segment. A respective software digital filter 417 
receives complex time domain output from the second FIFO DMA buffer 425' and 
performs the operations of tune and filtering in the complex frequency domain. 
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Therefore, each software synthesis receiver employs a respective FFT device 450, 
a respective filter Hm(f) 456, a respective step tune selector 452 employing offset indices 
for multiplier device 451 that multiplies filter weights of filter Hm(f) 456 and spectral 
domain data from FF T device 450 and outputs filtered and tuned data to IFFT device 453 
5 returning to the time domain. Data is output to respective I&Q phase rotation device 454 

to compensate for step tune phasing and connected to FIFO DMA buffer 455 for output. 
Respective two stage synthesis filter 416 and plurality of synthesis filter receivers 417 
are under control of the synthesis receiver control and demodulator processor 435. 

The synthesis receiver controller 435 activates, and controls two-stage synthesis 
10 filter 416 and activates, controls, and receives the output from one or more of the 

plurality of software synthesis receiver filters 417, performing signal demodulation in 
accordance with control instructions received from the real time controller 413. 

Data Buffering Device 

The demand memory access buffer 230, 330, 425, 425', 434 and the demand 

1 5 memory access buffer 412 are implemented as data buffering devices designed to receive 

synchronous or asynchronous data input such as from a previous module high-speed data 
router device. For efficiency, some modules of the system can be configured to share a 
common data buffering device when the data source for each of that module's 
components are the same. In general, these data buffering devices comprise a memory 

20 storage means such as computer processor RAM, cache, or other dedicated computer 

memory. These data buffer device are formed to index memory by First In First Out 
(FIFO) means, First In Tap Out (FITO) means, or data decimation means to provide 
asynchronous data output according to demand of the respective module digital 
processing device. Data may be input to and/or output from the buffering device as 

25 individual data or data blocks of various size. The same data may be asynchronously 

and/or time delay outputted to respective receiving modules. The output of the buffering 
device may represent time-domain data samples, time-decimated data samples, 
frequency-domain spectra data, or parameters domain data (for example angle-of-arrival 
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data per hyperchannel). The data buffering device may also provide selectable data 
indexing for data decimation wherein output data represents every mth data element. 

High-speed Data Router Device 

The high-speed data router 70 (not shown) is comprised in part as disclosed above 

5 with one to one, one to many, and many to one data router ports (DRP), exemplified by 

DRP 270 connected to plurality DRP 307, plurality DRP 370 connected to 407, as 
disclosed above, are used to distribute the output of the respective digital processor 
devices in which they are incorporated. As disclosed for the preferred described above, 
the data router port 270 is incorporated in the sub-band module 20 to output decimated 

10 wideband receiver data, while the data router port 307 is used to input respective data to 

the channelizer module 30. In similar fashion, the plurality of data router ports 370 of 
channelizer module 30 output data to plurality of data router ports 407 of signal processor 
module 40. In the embodiment of the present invention high-speed data output may be 
directed to another digital processing module, to a logging or archiving device, to a 

15 recording device, to a display device, or to other unspecified function device or module. 

Recursive routing of output data to input of same respective module may be performed. 
High-speed data routing is defined as the device capability to transfer data from one 
module to another without sustained and increasing backlog of data waiting to be 
transferred. 

20 High-speed data routing occurs at transfer rates substantially higher than the 

collective data rates or subsequent module output data rates. In one embodiment of the 
present invention, a commercial high-speed RACEway Interlink may be used to 
interconnect selected digital signal processing modules. The transfer route of data is 
configurable by the system controller 80 that is connected to the high-speed data router 

25 controller 70 (not shown) and to respective DRPs by the system control bus 82 (e.g. 

VMEbus). 

Transferred data may have various forms such as complex value in-phase and 
quadrature (I&Q) time domain data, complex value spectral data, baseband or audio 
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sample data, demodulated signal data, binary or symbol text data, or subsequent module 
notification data, such as signal activity frequency and bandwidth, from the detection 
processors 4 1 0 5 4 1 0' 9 4 1 (Teamed to respective demodulation/recognition processors 415, 
415^,415^ 

5 System Controller 

The system controller 80 is implemented in at least one embodiment using a VME 
control computer using such real time operating systems as Vx Works. The primary 
function of the system controller 80 is to coordinate module functions. The control 
function provides the overall set-up, coordination, management, and monitoring of the 

10 other system functions. Control sets up the individual functions by initializing the 

processes, downloading appropriate software and configuring the functions with 
necessary parameters. These parameters are typically maintained in a database and 
include data such as bandwidths, resolutions, thresholds, function sequencing, tasking, 
and processing selection. Control provides the system coordination function by 

1 5 implementing the defined processing sequences. This management allows for flexibility 

to soft-route data flows between functions. Management allows for the load balancing, 
resetting, and communications between functions. By monitor the individual functions 
control performs health check functionality. 

The real time controller 413, 413 5 , 413" is implemented using embedded 
20 processors such as the Motorola G3 processor and respective native operating system to 

control on-board module processes such as detection. The real time controller 413, 41 3\ 
413" is responsible for spawning software modules such as signal recognition and 
demodulation and direction finding. 

Signal Activity Detection 
25 One variation to the detection processor 414, as shown in Figure 4, involves the 

incorporation of a signal activity detection, the signal parameter encoder 443 that is 
connected to receive the moving average data from the buffers AVI - AV3 of the 
detection processor 414. As a further variation to the present invention, the use of the 
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signal parameter encoder 443 includes the generation of a long-term moving average 
LTA via a buffer AV3, along with the STA, MTA and DLTA averages. The signal 
parameter encoder 443 incorporates a buffers AV1-AV3 that receives the STA, MTA, 
and LTA or DLTA hyperchannel averaged data. A signal activity comparator (not 
5 shown) within signal parameter encoder 443 compares the STA, MTA and LTA or 

DLTA of each hyperchannel for a number of averages and compares to a predetermined 
threshold value to determine whether signal activity is present. The onset of a signal 
stronger than the ambient noise creates a temporary discrepancy between said STA and 
the LTA or DLTA in one or more adjacent hyperchannels encompassing the bandwidth 

10 of the signal activity. The signal activity detection (not shown) within parameter encoder 

443 characterizes time and amplitude of adjacent hyperchannels in the region surrounding 
a threshold excess. In operation, to mitigate the effects of pre-existing signals on the 
characterization of the new signal, the signal activity detection parameter encoder 443 
examines the ratio of MTA/LTA or MTA/DLTA. In further variations of the invention, 

1 5 the detection processor 414 may be devised to generate the STA, MTA, LTA and DLTA, 

not just the three moving averages used in the first embodiment. The STA, MTA and 
LTA or DLTA may in turn be used in other detection algorithms known to those skilled 
in the art to determine regions of hyperchannel signal activity. 

The signal activity detection parameter encoder device 443 alerts real-time 
20 controller 41 3 to spawn subsequent components or modules connected thereto to respond 

to detected and refine hyperchannel signal activity. Time, frequency and the number of 
hyperchannels (bandwidth) of signal activity exceeding predetermined thresholds may 
be distributed via high-speed routing devices (not shown) to such subsequent spawned 
processing modules for, as an example, synthesis, recognition, demodulation, and angle 
25 of arrival functions. 

Signal Recognition 

In the current embodiment of the invention, as shown in Figure 3 and Figure 4, 
channel processor 410 may be implemented using a plurality of dynamically allocated 
signal recognition devices 415, each respective signal recognition device comprising two- 
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stage synthesis filter 416 and plurality of synthesis filter receivers 417. Two-stage 
synthesis filter 416 receives complex spectral data from FITO data buffering device 412 
with a finite but variable delay tap of data. The FITO buffer 412 provides selected 
(decimated) complex hyperchannel data upon direction from a signal detection means, 
5 such as the signal detection device 414. Respective ones of the FITO buffer 412 may be 

configured to provide asynchronous data output to respective synthesis filter devices. It 
is to be understood that a plurality of synthesis filter devices may be spawned and 
controlled by a signal detection device to perform signal filtering on detected 
hyperchannel spectral activity using the same FITO memory or buffer data source. 

10 A synthesis filter device as would be applicable to the invention is based on the 

conventionally-known overlap-save method of linear FIR filtering. The basic operations 
in the overlap-save method are: (1) compute an N-point FFT of a complex time domain 
signal; (2) multiply it by the frequency domain representation of a filter; (3) compute an 
N-point inverse FFT of the product. These operations result in a filtered time domain 

15 version of the original signal. 

In the present invention, a number of significant changes to the basic overlap-save 
method are made. First, the data representation at the input of the two-stage synthesis 
filter 416 is a subset of FFT points of hyperchannel data originating from the N-point 
FFT device 310. Second, the multiply filter 420, 420' is comprised of real weight filter 

20 coefficients selected such that the filter product at multiply device 421 has 90 dB of 

attenuation at the band edge of the subset of FFT points and allows high signal fidelity 
to be achieved in the presence of strong (+80 dB) interference within the input 
bandwidth. Third, the multiply at 421 is offset in indices by step tune device 422. The 
shift caused by 422 effectively tunes the filter. The tuning step not only offsets the FFT 

25 data, but introduces a phase rotation. Respective tuned and filtered data is time domain 

transformed by inverse FFT (IFFT) 423. The resultant complex time domain data is 
contra-rotated by rotation device 424 to compensate for phase angle introduced by step 
tune device 422. 
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Another innovation of the present invention is the application of a second stage 
overlap-save operation to compute an N'-point FFT of a complex time domain signal; 
multiply it by the frequency domain representation of a filter; and return data to 
respective time domain data. FFT device 426 transforms time domain data received 
5 from FIFO 425 into complex spectral data that is multiplied by H2(f) filter weights 420' 

in multiplier device 421'. Data from multiplier device 421' is transformed to time 
domain by IFFT device 423' and connected to FIFO DMA buffer 425' for asynchronous 
output. To get the overall bandwidth reduction typically needed in the HF radio 
spectrum (e.g. from 2.7 MHZ to 3 kHz), the present invention used two stages of overlap- 
10 save filtering as described above. 

The invention provides significant processing efficiencies by the repeated 
*|j application of this synthesis filtering in stages. For a 1000 to 1 bandwidth reduction, the 

i|? present embodiment of the invention uses two stages of overlap-save synthesis filtering. 

% J In further embodiments or variations of the invention, multiple stages of synthesis 

15 filtering can be provided to attain the desired bandwidth reduction. Rotation of the 

» overlapped data also causes a phase discontinuity at each output block boundary that is 

H i corrected before subsequent processing is performed. The process of cascading synthesis 

J** filters is illustrated in Figure 5, employing two-stage synthesis filters 416 with respective 

time domain output received by synthesis filter 417. In like manner, time domain output 
20 of synthesis filter 417 is received by cascaded synthesis filter 417'. The cascade may 

continue for as many stages as necessary for data filtering, as exemplified by synthesis 
filter 417". At each cascade, complex time data may be output to digital recording 
means?. 

As disclosed above, the main embodiment of synthesis filtering done by the 
25 present invention performs two stages of synthesis filtering to convert time overlapped 

hyperchannels to time-domain data that more closely match the spectra of the detected 
signal activity of interest. The following illustrates the performance of the signal 
demodulation/recognition processor of the present invention in connection with the use 
of a signal detection device. 
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The signal detection device 414 alerts real-time computer 413 of signal activity 
in 1 1 adjacent hyperchannels out of an 8 192 -point FFT from hyperchannelizer FFT 310. 
Real-time computer 413 spawns demodulation/recognition processor 415 comprised of 
synthesis filter components. The signal synthesis components of the 
5 demodulation/recognition processor 415 receive 128-points of spectral signal activity 

encompassing the 1 1 adjacent hyperchannels of signal activity from FITO delay buffer 
412 corresponding to a narrowband signal detection event. First stage of two-stage 
synthesis filter component 416 receives the delayed FITO buffer data of surrounding 
signal activity of 128 hyperchannels. At first stage of the synthesis filtering 416, step 

10 tune selector 422 changes indices to represent center 1 1 hyperchannels of signal activity 

as zero Hz and then performs a multiply by filter weights 420 representing a passband 
of 5 kHz and a stopband (-90 dB) of 42 kHz. The synthesis filtering then performs a 128 
point Inverse FFT (IFFT). The output of this first stage of synthesis filtering decimates 
the data rate by 64 (8192 point FFT / 128 point IFFT). A second stage synthesis filtering 

15 performs 1024 point FFT overlapped by 25% at FFT device 426. The center 256 

frequency channels are multiplied by filter weights 420' representing a passband of 5 kHz 
and a stopband (-90dB) of 6 kHz, followed by a 256 point IFFT. The output of this 
second stage of synthesis filtering decimates the data rate by 4 (1024 point FFT / 256 
point IFFT) to a final complex rate of 10546.875. The synthesis filtering components 

20 then output time-domain data to a plurality of synthesis filter receivers 41 7 for refined 

signal data to the demodulation processor 435. 

The filters described in the synthesis processing are typically pre-computed and 
stored in a library. Different bandwidth processing may occur on the same initial FFT 
data simultaneously. For instance, in addition to the 5 kHz bandwidth filter described, 
25 a simultaneous two-stage 50 kHz synthesis filter may be instantiated. The first stage of 

this 50 kHz synthesis filter processing would use a passband filter of 50 kHz and a 
stopband of 84 kHz, followed by a 256 IFFT. The second stage synthesis filter would 
then consist of a 512 point FFT overlapped by 12.5%. The second stage filter would 
have a passband of 50 kHz and a 90 dB stopband of 56 kHz and final 512 point IFFT. 
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Synthesis Receiver Control and Demodulation Processor 

The synthesis receiver control and demodulation processor 435 receives a 
decimated and synthesis filtered data stream from the synthesis filter components 41 7 of 
the demodulation/ recognition processor 415. 

5 The processor, as shown in Figure 4, incorporates a signal recognition data buffer 

425' that provides the data stream to a plurality of signal synthesis filter receiver devices 
417 that performs filtering in a manner similar to a synthesis filter module. To continue 
illustration of the current embodiment responding to a detection event of 11 
hyperchannels, the signal recognition device receives the 10546.875 Hz sample rate data 
10 center tuned to the frequency of signal activity. The signal synthesis filter receiver 

performs a 512 point, 25% overlapped FFT and convolves with the three non-zero 
coefficients of the DFT of a Hann window. A resolution of about 40 Hz is thereby 
achieved. The original hyperchannel signal activity is thus resolved into spectral 
components. 

15 In addition, the signal synthesis filter control and demodulation processor 435 

makes a run-time decision about the need for additional filtering. As a typical example, 
given two spectral components indicating a Frequency Shift Key (FSK) signal, the 
synthesis filter control and demodulation processor applies three simultaneous single 
stage synthesis filters 417 one on each individual spectral components and one on the 

20 total apparent signal of interest bandwidth. The criteria for the number, placement, and 

span of synthesis filters in the hyperchannelized data is maintained in a controlling 
database (not shown) and initiated by detection processor 414 through real-time computer 
controller 413 to synthesis receiver control and demodulation processor 435. 

The synthesis receiver control and demodulation processor 435 applies multiple 
25 and simultaneous signal demodulations to the time domain data output from plurality of 

synthesis filter receivers 417. The algorithm performs multiple and simultaneous 
demodulations based on different signal modulation types and keying speed. Synthesis 
receiver control and demodulation processor 435 performs single-channel amplitude, 
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frequency, or phase demodulation by generic algorithms in conjunction with a database 
of signal demodulation criteria. The criteria may be varied to control probability of 
correct signal identification and probability of false alarm. 

In another variation of the invention, the synthesis receiver control and 
5 demodulation processor 435 makes a run-time decision that a signal is compound 

modulated, such as by frequency modulated frequency division multiplexing (FM/FDM) 
of narrowband telephony circuits. To accomplish this, the signal recognition processor 
performs FM demodulation on the output of the first synthesis filtering. First synthesis 
filter bandwidth is examined by synthesis receiver control and demodulation processor 

10 435 and determines that, for example, 12 or more narrowband frequency division 

multiplexed signal channels could be further demodulated, synthesis filter control and 
demodulation processor 435 then subsequently applies 12 simultaneous single stage 
synthesis filters 417, one per channel. The digital output from respective synthesis filters 
are received by synthesis filter control and demodulation processor 435 and applies 

15 generic algorithms in conjunction with a database of signal demodulation criteria and 

performs single-channel amplitude, frequency, or phase modulation. The criteria may 
be varied to control probability of correct signal identification and probability of false 
alarm. 

N-Channel Direction Finding 

20 It is to be understood that the present invention may employ wideband data from 

a plurality of antennas or sensor and receive multiple respective wideband data. 
Wideband data may be received from a set of N spatially separated or phase orthogonal 
antennas or sensors, and as is common for such multi-antenna devices, the respective 
multiple of wideband receivers are coherently phase locked and ADCs are time 

25 synchronized. 

In accordance with a second embodiment of the invention, as illustrated in Figure 
6 for three spatially separated or phase orthogonal antennas or sensors 5, 5\ 5", 
connected to a respective plurality of multiple wideband receivers 110, 110', 110" that 
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have coherently phase locked oscillators (not shown) and output coherently phased signal 
data to sub-band decimation module 20 that receives respective time domain analog 
signals and converts respective sensor signals into complex digital data outputs from a 
plurality of multiple of wideband ADCs 210, 210', 210". In this embodiment of the 
invention, six decimated complex digital sub-band data streams are output from each 
respective antenna or sensor source wideband ADC. Figure 6 illustrates this respective 
plurality of sub-band data streams exemplified by two sub-band decimation devices 220 
with I&Q time domain output of respective sensor 5. In similar fashion, two sub-band 
decimation 220' and 220" outputs of respective sensor 5' and 5" are illustrated. 

The plurality of multiple wideband decimation outputs 270, 270', 270" are 
distributed to a respective plurality of FFT hyperchannelization modules 3 10, 3 10', 310", 
connected to respective FIFO DMA buffer devices 330, 330', 330" and output distributed 
by high-speed data router ports 370, 370', 370". 

The object of the FFT devices 310, 310', 310" is to channelize respective N- 
sensor time domain data into complex spectral components that are much narrower than 
any signal of interest bandwidth. FFT data from respective antenna or sensor 
multichannels are time synchronized. 

In this second embodiment of the invention, a first channel of FFT hyper 
channelized data 370 provides input to the spectral detection circuit 414 of first channel 
processor 410 by means of data router port 407. Upon detection of signal activity, the 
spectral detection circuit spawns the operation of a demodulation recognition processor 
415 comprising synthesis filter 416, a plurality of simultaneous synthesis filter receivers 
417 whose respective time domain signal output is demodulated by synthesis receiver 
control and demodulation processor 435 similar to those used in the first embodiment of 
the invention. With N multichannel antenna or sensor source FFT hyperchannel data, the 
spectral detection circuit 414 of first channel processor 410 may also spawn a direction 
finding processes within direction finding module 60. 
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One implementation of the direction finding module 60 involves the use of the 
N-Channel Direction Finding technique. Referring to Figure 6 and 8a, DF processor 
module 60 comprises a plurality of input data router ports 607 and 607' , receiving upper 
and lower half hyperchannelized data from respective decimated sub-band of N-sensor 
5 channel processors. Decimated sub-band hyperchannelized data of N-sensor channel 

processors from respective data router ports 607 and 607' are collectively connected to 
respective DF Processor 610. DF Processor 610 comprises respective upper and lower 
hyperchannel data multiplexers (MUX) 611,611' connected to respective plurality FITO 
delay memory buffers 612, 612'. In response to detected signal activity by detection 

10 processor 414, respective N-sensor hyperchannel data of FITO delay memory buffers is 

output to respective spawned N-sensor band DF two-stage synthesis filters 614, 614'. 
Filtered N-sensor upper and lower half hyperchannelized data is transformed by 
respective DF two-stage synthesis filters 614, 614' to time-domain data and output to 
respective spawned upper and lower half hyperchannel DF calculators 616, 616' to 

15 perform angle of arrival calculation. DF calculators 616, 616' calculate angle of arrival 

of signal using N-sensor collective time domain data by phase comparison algorithms 
known to those skilled in the art. Calculated angle of arrival parameter is output to data 
bus 83. 

N-Channel Multiplex 

20 N-sensor FFT channel processors 310, 310', 310" is distributed by high-speed 

data router and FIFO buffers to upper half and lower half N-sensor data multiplex devices 
611,611' of DF processor 610. It is to be understood that the hyperchannel complex 
spectral data corresponding to the upper half spectral data of respective sub-band 
decimation FFT channel processor of each respective sensor is collectively received by 

25 upper half multiplexer 611. Respective hyperchannel complex spectral data 

corresponding to the lower half spectral data of respective sub-band decimation FFT 
channel processor of each respective sensor is collectively received by lower half 
multiplex 611'. 
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As an example, a 4-sensor system comprises four respective wideband receivers 
110 output to four respective ADCs 210. Each respective ADC 310 outputs to six sub- 
band digital down converters 220. Time domain data from plurality of 24 digital down 
converters 220 is routed by high-speed data router to respective plurality of 24 FFT 
channelizers 3 1 0 transforming data into hyperchannel complex spectral data that is output 
and distributed by high-speed data router to respective six upper half multiplexers 611 
and to respective six lower half multiplexers 611', wherein each multiplexer receives 
respective sub-band decimated hyperchannel data from each of the four sensors. 
Referring to Figure 10a, respective multiplexer device receives hyperchannel complex 
spectral data from N-sensor FFT channelizers and outputs N-sensor data streams to 
respective FIFO delay memory data buffers. 

DF Synthesis Filter 

The multichannel DF synthesis filter 614 performs N-sensor simultaneous 
synthesis filtering on a common respective subset of hyperchannels that are centered on 
the hyperchannel of first channel detected signal activity. In the present embodiment of 
the invention, two-stage synthesis filtering is used in a manner similar to the first channel 
synthesis filter. Respective synthesis filters perform identical bandwidth filtering and 
provide N multichannel filtered time domain outputs of detected signal activity. 
Referring to Figure 9a, N-sensor hyperspectral data output from plurality of FITO delay 
memory buffers 612a, 612b, ... 612n is received by respective first stage multiple 
synthesis filter multiplier devices 621a, 621b, ... 621n of spawned DF two-stage synthesis 
filter 614. Respective multipliers multiply hyperchannelized spectral data with common 
frequency-domain filter coefficient Ho(f) table 620 using common step-tune selector 
offset indices 622. In manner similar to first stage synthesis filter 416, frequency filtered 
hyperchannel data from multiplier devices 621a, 621b, ... 62 In are output to respective 
inverse FFT devices 623a, 623b, ... 623n connected to respective contra-rotating phase 
rotation devices 624a, 624b, ... 624n connected to respective FIFO DMA buffers 625a, 
625b, ... 62 5 n for asynchronous output respective second stage filters. FIFO DMA 
buffers 625a, 625b, ... 625n output time-domain data to respective FFT devices 626a, 
626b, ... 626n that transform N-sensor data back to complex spectral data and output to 
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respective second stage multiplier devices 627a, 627b, ... 627n multiplying respective data 
by common second stage filter coefficients Ho'(f) 620'. In manner similar to first stage 
of multiple N-sensor synthesis filters, multiplier devices 627a, 627b, ...627n output is 
received by respective inverse FFT devices 623a', 623b', ... 623n' that transform complex 
spectral data back to time-domain data received by respective FIFO DMA buffers 625a', 
625b' , ... 625n 5 for asynchronous output of N-sensor tuned and filtered time-domain data. 

DF Calculator 

The DF angle of arrival calculation device 616 receives sensor time domain data 
for calculating an angle of arrival (AOA). As will be evident to those skilled in the art 
of direction finding, any of a number of commonly known phase and amplitude N sensor 
comparison algorithms may be employed by the DF angle of arrival calculation device 
616 of the present invention. For example, in this second embodiment of the invention, 
an angle of arrival calculation device 616 would be able to accommodate 2-channel 
commutated direction finding calculation, 3 -channel comparative direction finding 
calculation, as well as any of a set of commonly known N multichannel comparison 
algorithms. 

3 -Channel Comparative DF Module 

A further implementation of the direction finding embodiment uses the 3 -Channel 
Comparative DF technique. In this implementation, referring again to Figure 6, three 
sensors are employed: reference antenna 5, sine-phase antenna 5\ and cosine-phase 
antenna 5" provide the RF or acoustic source. A first channel reference source is 
received by the first wideband receiver 1 1 0 and ADC 210. The second and third channel 
sources are received, respective, by the second and third wideband receivers 1 10', 110" 
and ADCs 210', 210". The first, second, and third local oscillators of wideband receivers 
1 10, 1 10', 1 10" are coherently phase locked and the corresponding ADCs 210, 210, 210" 
are time synchronized. 

First, second, and third channel digital data are received by the respective 
wideband decimation modules 220, 220', 220" and respective FFT hyperchannelization 
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modules 310, 310% 310". First channel hyperchannel data is received by first channel 
signal detection processor 414 of channel processor 410. The spectral detection 
processor 414 spawns the operation of the synthesis filtering module 415, comprising 
two-stage synthesis filter 416, a small plurality of synthesis filter receivers 416 and 
synthesis receiver control and demodulation processor 435. Concurrently, respective 
first, second and third channel hyperchannel data are received by 3 -channel comparative 
DF processor 610. 

Referring to Figure 8b, upper and lower spectral half data of respective first, 
second, and third FFT hyperchannel data are received by respective multiplexer devices 
61 1 and 611 \ The three channels of data pass through the respective multiplexer as 
shown in Figure 10b, representing a 3-sensor instantiation of an N-sensor data flow. 
Again referring to Figure 8b, upper and lower half spectral data from multiplexer devices 
61 1, 611 5 are output to respective two-stage DF synthesis filters 614, 614' that output 
three streams of identically filtered and transformed time domain data to respective DF 
calculators 616, 616' 

The multichannel two-stage DF synthesis filters 614, 614' performs three 
simultaneous synthesis filterings on respective reference, sine-phase and cosine-phase 
respective subsets of hyperchannels that are center tuned on the hyperchannel of first 
channel detected signal activity. In this 3 -channel comparative DF implementation, two- 
stage synthesis filtering is used in a manner similar to first channel recognition synthesis 
filter. Respective synthesis filters perform identical bandwidth filtering and provide three 
filtered time domain output of detected signal activity for output to DF calculator 616, 
616'. 

In another variation of this implementation, additional synthesis filters may be 
incorporated to perform fine grain spectral filtering on the three comparative channel data 
based on the first channel synthesis filters of the signal recognition module 414. 
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An angle of arrival calculation device 6 1 6, 6 1 6* receives reference, sine-phase and 
cosine-phase time domain and filtered data of comparative sensor hyperchannels. As one 
of skill in the art would understand, amplitude and phase angle correlation, Watson- Watt 
amplitude, or other angle of arrival algorithms may be used to calculate the angle of 
5 arrival. Angle of arrival data is output to data bus 83. 

2-Channel Commutated Direction Finding Module 

With respect to the fourth embodiment of the invention, it is to be understood that 
the present invention is also capable of implementing specific multichannel direction 
finding (DF) algorithms. One such implementation involves the 2-Channel Commutated 
10 DF technique. 

As shown in Figure 7, reference antenna 5, sine-phase antenna 5', and cosine- 
phase antenna 5" provide RP or acoustic data sources to a dual wideband receiver 111. 
A first channel reference source 5 is received by a first wideband receiver of dual 
receivers 1 1 1 and a first ADC 210. A second channel cyclically commutates between 
15 reference 5 , sine-phase 5' and cosine-phase 5" sources as input to a second wideband 

receiver of dual receivers 111 and a second ADC 210'. The first and second wideband 
receivers 1 1 1 are coherently phase locked, and the first and second ADCs 210, 210' are 
time synchronized. 

First and second channel digital data output from the first and second respective 
20 ADCs 210, 210' are received by respective plurality of wideband decimation devices 

220, 220' and respective FFT hyperchannelization processors 310, 310\ The first 
channel of hyperchannel data is received by respective plurality first channel of channel 
processors 410 comprising respective plurality of first channel spectral detection 
processor 414. Respective first channel spectral detection processors 414 spawn the 
25 operation of first channel demodulation synthesis filter processor 415 comprising two- 

stage synthesis filtering device 416, a plurality of synthesis filter receivers 417, and 
synthesis filter control and demodulation processor 435 to detect, demodulate and 
identify a signal of interest. Concurrently, respective first channel and second 
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commutated channels of hyperchannel data are received by respective plurality of 2- 
channel comparative DF processors 610. 

Referring to Figures 7 and 8c, respective data port router 607 and 607' of DF 
processor module 60 receives respective upper and lower spectral half of sub-band of 
first and second FFT hyperchannel data and outputs to respective DF processor 610. DF 
processor 610 comprises multiplexers 61 1 and 61 1 5 that receive respective upper and 
lower half data and output data to respective three channel upper and lower half FITO 
delay memory buffers 612, 612' that output hyperchannelized data from reference 
channel, sine-phase channel, and cosine-phase channel to respective upper and lower half 
DF synthesis filters 614, 614' connected to respective DF calculators 616, 616' for AOA 
output to data bus 83. As one of skill in the art would understand, amplitude and phase 
angle correlation, Watson- Watt amplitude, or other angle of arrival algorithms may be 
implemented by DF calculator to calculate the angle of arrival. 

2-Channel DF Multiplexer 

Referring to Figure 10c, it can be seen that multiplexer device for receiving 2- 
channel reference and commutated hyperchannel data operates differently from 
multiplexers employed by N-channel or 3-channel comparative DF operation. Where as 
N-channel and 3-channel data is collectively received by multiplexer 61 1 and directly 
output to N plurality or 3 plurality FITO delay memory buffers as illustrated in Figures 
10a and 10b, 2-Channel DF multiplexer 611 operation as illustrated in Figure 10c 
decommutates second channel FFT hyperchannel commutated data. In the embodiment 
of the present invention incorporating 2-Channel DF multiplexer 611 performs two 
operations on second channel FFT hyperchannel commutated data: (1) time commutated 
second channel data is separated into three channels representing reference source, sine- 
phase source, and cosine-phase source hyperchannel data; (2) FFT hyperchannel 
commutated data calibration is performed. 

2-Channel DF multiplexer performs calibration by operation of comparing 
respective first channel reference data sample and second channel reference data samples 
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representing synchronous FFT hyperchannel data Since identical synchronous reference 
source 5 is origin of both first and second channel hyperchannel are compared by 
multiplexer operation, difference between first and second channel reference data is 
attributed to difference in analog filters (not shown) and dual wideband receivers 111. 
Therefore same difference is applied to following second channel sine-phase and cosine- 
phase FFT hyperchannel samples. A calibration operation in 2-Channel DF multiplexer 
611 performs adjustment on respective commutated channel data in such a manner that 
all hyperchannels of the second channel are adjusted an amount to ensure the reference 
data of the second channel equals the reference data of the first channel. 2-Channel 
multiplexer 61 1 therefore applies first to second channel reference difference to following 
second channel reference, sine-phase channel, and cosine-phase channel FFT 
hyperchannel data. 

2-Channel DF multiplexers 611,611' outputs corrected reference channel, sine- 
phase channel, and cosine-phase FFT hyperchannel data to respective three FITO delay 
memory buffers 612, 612'. 

2-Channel DF Synthesis Filters 

The DF synthesis filtering device 614, 614' performs three synthesis filters on the 
plesio- synchronous second channel reference, sine-phase and cosine-phase respective 
subset of hyperchannels output from respective 2-Channel FITO buffers 612, 612' 
received from 2-Channel DF multiplexers 611,611'. In this 2-channel commutated DF 
embodiment, two-stage synthesis filtering is used in a manner similar to the first channel 
synthesis and provide three filtered time domain outputs of detected signal activity. 

In yet another implementation of the direction finding embodiment, additional 
cyclic synthesis filters can be added to perform fine grain spectral filtering on the second 
channel data based on the first channel synthesis filtering of the signal demodulation and 
recognition processor 415. Specifically, multichannel data is filtered in like manner of 
first stage of two-stage DF synthesis filter 614 comprising plurality of channel multipliers 
connected to common shape filter and common step tune selector multiplier indices 
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offset. Multichannel multiplier product is output to respective multichannel inverse FFT 
connected to respective multichannel phase rotators connected to respective multichannel 
FIFO DMA buffers for output. In this preferred embodiment, output from plurality of 
fine-grain synthesis filters for each source channel are received by plurality of DF 
calculators. AOA for each respective fine-grain syntheses filter is computed according 
to comparative DF algorithm or other algorithms known to those skilled in the art. DF 
calculator outputs plurality of AOAs from respective plurality of fine grain synthesis 
filters to data bus 83. Subsequent data analysis processes may use plurality of AOAs 
from signal activity as individual representations of direction or may be averaged to 
provide statistical representation of direction of arrival as performed by well understood 
mathematical averaging operations. 

Other Direction Finding Embodiments 

In the above-described implementations of direction finding embodiment of the 
invention, a received subset of delayed FFT hyperchannel data from each antenna or data 
source is selected by a first channel spectral detection processor and used to spawn 2- 
channel commutated DF, 3-channel comparative DF or N-Channel DF synthesis filter 
and DF calculation processor for angle of arrival calculation. 

However, in other potential DF embodiments of the present invention, the 
selected subset of FFT hyperchannel data may be known a priori for fixed channel signal 
activity and may be immediately selected for DF synthesis filtering and DF calculation, 
and need not be dependent upon first channel spectral detection or signal recognition 
process direction. Alternatively, the complete set of FFT hyperchannel data may be 
accepted from N multichannel sources. Here, synthesis filtering is not required, and a N- 
multichannel comparative or 2-Channel commutated DF algorithm is performed on each 
hyperchannel. 

It is to be understood that, in the DF embodiments of the present invention, 
performing angle of arrival determination on every hyperchannel, a direction finding 
module may precede and replace a signal recognition and demodulation module 
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comprising signal detection processor. An alternative signal activity detection processor 
would be used to compare multiple adjacent hyperchannels for amplitude and angle of 
arrival consistency. Amplitude and angle of arrival consistency criteria for various 
signals of interest are simultaneously applied by respective DF consistency signal activity 
5 detection processor. As with the spectral detection processor, a DF consistency signal 

activity detection processor spawns the operation of module processes of synthesis 
filtering and signal recognition. In this embodiment of the invention, data 83 is replaced 
by high-speed data routers to output signal activity hyperchannel center (frequency), 
number of adjacent hyperchannels (bandwidth), and time of hyperchannel signal activity 
10 to respectively spawned digital signal processing modules. 

For summary of the above embodiments, Figure 11a depicts arrangement of 
receiver and processing modules for single sensor source sub-band decimation, signal 
hyperchannelization, signal recognition and demodulation; Figure lib depicts 
arrangement of N sensors 5 ... 5", connected to corresponding N receiver modules, first 

15 channel processing modules of sub-band decimation, signal hyperchannelization, signal 

recognition and demodulation, N-channels of sub-band decimation, signal 
hyperchannelization and collective N-channel module for direction finding (including N- 
sensor and 3 -sensor comparative data for DF calculation); Figure 1 lc depicts invention 
embodiment of 3 sensor sources with commutated second channel, arrangement of first 

20 channel receivers and respective first channel processing modules of sub-band 

decimation, signal hyperchannelization, signal recognition and demodulation, and second 
channel receiver and second channel processing modules of sub-band decimation, signal 
hyperchannelization, and collective second channel module for direction finding (using 
calibrated 3-sensor commutated data for DF calculation); Figure 1 Id depicts N-channel 

25 invention embodiment wherein DF processor precedes signal detection and recognition. 

DF information from collective N-channel DF processor is output by high speed data 
router to first channel signal processor, performing signal activity detection and spawning 
of signal recognition and demodulation processors based upon consistent angle of arrival 
direction over desired signal bandwidth. 
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Beamforming 

Yet another embodiment of the invention is to perform beamforming to enhance 
initial signal detection. Referring to Figure lie, processing modules are again re- 
arranged. In this embodiment, beamformer module 690 provides a means for spatial 
5 enhancement of N-sensor data prior to signal detection modules 40. Beamforming sensor 

data may be performed on a wideband or narrowband plurality of hyperchannel data. As 
illustrated in Figure 8d and 9b, one or more plurality of beamforming processors 617 use 
beamform weighting 630. Beamformer 617 receives N-sensor hyperchannel data from 
respective FITO delay memory buffers 612a to 61 2n and performs filtering and tuning 

10 connected to respective inverse FFT 623a to 623n and phase rotation compensators 624a 

to 624n that output time-domain data to respective FIFO data buffer 625a to 625 n. 
Respective N-sensor time-domain data from data buffers 625a to 625n are multiplied by 
N-sensor beamforming weights 630 by respective multiplier devices 631a to 63 In, 
summed by device 635, and output to FIFO buffer 636. Buffered beamformed summed 

1 5 output 636 is received and transformed to hyperchannel frequency domain by FFT device 

637 and output to high speed router device 70. 

Beamforming may be performed on any plurality of hyperchannels including the 
plurality of all hyperchannels representing the full bandwidth of the wideband front-end 
receivers and channelizer modules 30. Beamforming weights may be provided a priori 

20 as illustrated in Figure 9b. It is to be understood that amplitude and/or phase biases of 

the sensor and RF paths (i.e., sensor location, RF cables, RF distribution, filtering and 
digitization) are controlled, and characterized between all N sensor channels. A 
calibration table or other such normalization scheme characterizes constant biases. The 
calibration table is used to compensate for, and normalize each sensor RF path. The 

25 wideband data is digitized at a rate greater than Nyquist for the frequency of interest and 

routed to pre-processing functions. The digitization process must be done such that the 
phase of the digitization clocks remains constant between all sensor channels. 

Beamforming may be performed in conjunction with signal processor module 40, 
wherein the hyperchannels selected for beamforming represent the bandwidth of a 
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specific signal of interest. Beamforming may be performed prior to signal detection, 
from such common channel signals as cellular, wireless, and other personnel 
communication services (PCS) where signal activity is known a-priori in specified 
narrowband channels. 

To achieve the beamformed processing the N sensor hyperchannels are combined 
into a single path by means of a sum-of-weighted-taps known to one of ordinary skill in 
the art may be implemented either in the time domain as a sum-of-weighted-taps or in the 
equivalent FFT / hyper-channelization transfer function frequency domain processing. 

The N-sensor beamforming may be performed at various stages of processing (i.e. 
before decimation, before or within the hyper-channelization process, or before, within 
or after the matched bandwidth processing). 

The flexibility in order of processing stages provides for a means to balance and 
optimized process loading. The values of the sum-of-weighted-taps weights and taps 
provide the spatial filtering pattern. The values of the sum-of-weighted-taps weights and 
taps are applied across the band to which that spatial filtering is to be applied. The 
particular spatial filter may be applied to the complete spectra or one or many sub-bands. 
The band of application is known at the time of application. The band may be known by 
a priori information or as the result of another processing stage (i.e. Detection, 
Narrowband processing, etc.) The particular spatial filter being applied is also known at 
the time of application. The spatial filter may be known by sector of interest, or could be 
derived from another processing stage (i.e., DF or a previous Beamformed processing 
stage). The spatial filter and band of application may be altered or modified during 
processing or (because of data buffering) re-applied from a previous time. The path 
combination process function (at whatever stage included) provides the data to a router 
function for distribution to a follow-on processor stage. Each processing stage applied 
before path combination (whenever applied and including path combination) must take 
care to not introduce any biases between each of the N-sensor paths. This is done by 
taking care to performing identical processing to each respective sensor path and assuring 
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the identical processing functions are applied equivalently to the data samples for each 
path digitized at the same time. This requires that any transfer function applied (i.e. FIR 
filter, IIR filter, FFT, IFFT, decimating filter, hyper-channelization, cyclostationary 
process or the like) be started and operate on data samples generated from the same 
5 digitization clock cycle from each sensor path. For a typical embodiment, a pre- 

processing hyper-channelization function converts the time samples to a series of 
frequency spectra for each sensor path at a rate greater the Nyquist for the individual FFT 
bin resolution. The rate in which the FFT is greater than the Nyquist rate for that 
frequency bin is the overlap percentage. Typical overlap percentages for Beamedformed 
10 collection processing are between 6.25 to 87.5%. The size of the FFT bins is such that 

each is of finer resolution than an individual signal waveform. Typical FFT hyperchannel 
resolutions are 1 to 25% of signal waveforms. 

Adaptive Beamforming 

If the frequency and characteristics are not completely known a priori, a 

15 processing function may be applied to detect signal presents and estimate new or 

improved beamforming weight parameters including weights for either signal 
enhancement (beamforming) or weights for interference rejection (adaptive nulling). This 
process may also monitor control channels to determine similar information. A router 
function may then route the data in its entirety or as a subset to subsequent processing 

20 modules as illustrated in Figure lie, showing beamform hyperchannel data provided to 

signal processing module 40. 

Time Difference of Arrival (TDOA^ and Frequency Difference of Arrival (TDOA^ 

Time Difference of Arrival (TDOA) and Frequency Difference of Arrival 
(FDOA) signal processing provides a means for band limited hyperbolic direction of 
25 arrival and location processing that is match to the frequency of the signal energy or 

waveform. To perform TDOA processing, received signal or waveforms must be time- 
tagged in a consistent reference frame, allowing two or more plurality of systems to time 
compare signal data. In the embodiment of the present invention, analog to digital 
conversion in ADC 210 of sub-band module 20 applies time-tags to digitized signal or 
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waveform time domain data packets using an accurate time reference (not shown in 
Figure la or lb) such as a Global Positioning System (GPS) clock, cesium-beam clock, 
or other such accurate time keeping means. 

The time, group, amplitude and phase bias and/or delays of the RF paths (i.e. 
5 antenna location, RF cables, RF distribution, filtering and digitization), time 

synchronization (i.e. GPS receiver, time code generator, 1 pulse-per-second circuitry, 
digitization clock, and time reference,) and sample time association (i.e. time resolution, 
time representation, and time communication) are controlled and characterized for all 
channels. Ambiguity in the characteristics of the paths that are not known a priori may 
10 be calibrated. A calibration table characterizes all constant biases and said table is used 

to compensate for, and normalize each path. 

The wideband data is digitized by ADC 210 at a rate greater than Nyquist for the 
receiver bandwidth of interest and routed to a pre-processing function. For TDOA, two 
or more plurality of separate sensor systems as illustrated in Figure la comprising 
15 receiver module 10 connected to sub-band module 20 connected to channelizer module 

30 connected to signal processing module 40 are required to measure sensor data time of 
arrival from spatially diverse sensors. 

In another variation of the TDOA embodiment of the present invention, spatially 
diverse N-sensor data may be brought together in one N-channel system as illustrated in 

20 Figure lb, receiving N-sensor wideband data by N-plurality of receiver modules 10 to 

10" connected to respective sub-band modules 20 to 20" connected to respective 
channelizer modules 30 to 30" connected to respective signal processor modules 40 to 
40". First N-sensor channel signal processor module 40 performs signal detection and 
spawns identical signal demodulation and recognition processes in respective N-sensor 

25 signal processor modules 40 to 40". Common means of time of arrival measurement 

(e.g. time of first baud, linear regression of a set of bauds to estimate time of first baud, 
cross correlation of signal with a priori sequence, cross correlation of signal with first N- 
sensor channel signal data) from N-sensor channel signal processor modules 40 to 40" 
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are compared by first N-sensor channel signal processor module 40, calculating Time 
Difference of Arrival (TDOA) for signal activity. To achieve the TDOA processing, the 
channels are correlated and compared with data collected from the responses of the 
spatially diverse antennas to determine geolocation. 

The digitization process in each N-sensor channel must be done such that the 
phase of the digitization clocks remains constant between all channels. Time is associated 
with each digitization sample. This sample time associated remains linked to the data 
sample through data decimation, FFT hyperchannelization, FITO and FIFO buffering, 
filtering, tuning, and inverse FFT processing. 

Therefore, time of arrival measurement or correlation processing can be 
performed at various stages of processing (i.e. before decimation, before or after the 
hyper-channelization process, or before or after the matched bandwidth processing). The 
flexibility in order of processing stages provides for a means to balance and optimized 
process loading. The time of arrival or correlation process may be applied to the complete 
spectra or one or many sub-bands (possibly overlapping) or to synthesis filter bands. The 
band of application is known at the time of application. The band may be known by a 
priori information or as the result of another processing stage (i.e., detection, synthesis 
filtering). The integration time used for time of arrival measurement (e.g., linear 
regression) or correlation may be known in advance of processing or may be determined 
during processing (i.e., amount of convergence). The band of application and integration 
time of processing may be altered or modified during processing or (because of data 
buffering) and re-applied from a previous time. The TDOA process function (at whatever 
stage included) provides the data to a router function for distribution of time 
measurement data to a follow-on processor stage. The TDOA correlation may occur on 
bandlimited data samples or waveform features. Waveform feature correlation allows for 
efficiently passing less data between paths for those waveforms that have extractable 
features. Each processing stage applied before TDOA processing (whenever applied and 
including TDOA processing) must take care to not introduce any biases between each of 
the paths. This is done by taking care to apply identical processing to each and assuring 
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the identical processing is applied equivalently to the data samples for each path digitized 
at the same time. This requires that any transfer function applied (i.e. FIR filter, IIR filter, 
FFT, IFFT, decimating filter, hyper-channelization, cyclostationary process or the like) 
be started and operate on data samples generated from the same digitization clock cycle 
5 from each path and time association exists across each transfer function. 

For a typical embodiment, a pre-processing hyper-channelization function 
converts the time samples to a series of frequency spectra for each path while maintaining 
the sample time association across the transfer. A router function may be applied to this 
data to parse the data by FFT bins (frequency) or this function may be applied to the 
10 entire set of FFT bins. 

Additional stages of processing may be applied to optimized the band to be 
TDOA processed. A router function may be applied to this data. The router function 
(possible from spatially separate processes) brings the data streams from each of the 
path's data processing together for the TDOA time comparison or cross correlation 

15 process. The results of this process is a hyperbolic line of position and other parameters 

(i.e. quality, confidence, etc), or a time series of such based parameters. The resultant 
data may be combine in spatial or time-domain processing well known to those practiced 
in the art of TDOA geolocation to combine TDOA hyperbolic data into a geolocation and 
other parameters (i.e. quality, confidence, etc.) or a time series of such geolocation based 

20 parameters. 

Frequency Difference of Arrival may be combined with the TDOA processing by 
techniques known to one skilled in the art to contribute to the geolocation process. FDOA 
contributions can either replace a necessary TDOA path or supplement the set. The data 
may be routed to functions that repeat select processing techniques already applied. This 
25 may typically include additional hyper-channelization stages and additional spatial 

filtering to apply an even more closely matched filter to provide optimized frequency 
waveform processing. As is well known, wide bandwidth allows optimal time of arrival 
estimates and narrow bandwidth allows optimal frequency estimates. The embodiment 
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of the present invention can perform both optimizations by simultaneous processing paths 
with optimized hyperchannel bandwidths for both time correlation processes and 
synthesis filter frequency measurement processes. 

Although the present invention has been fully described in connection with the 
preferred embodiment thereof with reference to the accompanying drawings, it is to be 
noted that various changes and modifications will be apparent to those skilled in the art. 
Such changes and modifications are to be understood as included within the scope of the 
present invention as defined by the appended claims, unless they depart therefrom. In 
addition, the various embodiments and implementations of the present invention as 
disclosed and claimed above are suitable for use in many different industries and 
applications, including but not limited to the cellular, RADAR, medical diagnosis and 
data analysis, SONAR, acoustic, ultrasonic, video, broadcast, communications, 
telecommunications, geological oceanographic and natural science fields. 
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